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The Ziman formulation of electrical conductivity is tested in warm and hot dense matter us¬ 
ing the pseudo-atom molecular dynamics method. Several implementation options that have been 
widely used in the literature are systematically tested through a comparison to accurate but expen¬ 
sive Kohn-Sham density functional theory molecular dynamics (KS-DFT-MD) calculations. The 
comparison is made for several elements and mixtures and for a wide range of temperatures and 
densities, and reveals a preferred method that generally gives very good agreement with the KS- 
DFT-MD results, but at a fraction of the computational cost. 
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I. INTRODUCTION 

The calculation of electron transport coefficients in 
warm and hot dense matter has has been the subject 
of a considerable amount of recent theoretical interest 
[IH3 due in part to the importance of such quantities 
for inertial confinement fusion modeling [1, Q . A lack of 
direct experimental data for electrical and thermal con¬ 
ductivities also drives the need for accurate and reliable 
theoretical calculations. 

Kohn-Sham density functional theory molecular dy¬ 
namics^ has emerged as a powerful and accurate tool 
for simulating warm dense matter. By coupling it with 
the Kubo-Greenwood formulation [ol-fl^ it has also been 
used to generate conductivity data. Unfortunately, the 
method is extremely computationally expensive, and for 
this reason is typically limited to low temperatures (a few 
tens of eV). 

A much more computationally efficient alternative to 
QMD are average atom models. These come in many 
varieties [TsI - fT^ and are also density functional theory 
based models. The common theme to these models is 
that one tries to find the properties of one atom in the 
plasma. Average atom models have also been used to 
generate electronic transport coefhcients. Broadly speak¬ 
ing, there are to two main approaches: on the one hand 
are Kubo-Greenwood based calculations 0, iO-il that 
are qualitatively similar to the QMD calculations; and 
on the other hand are Ziman-type methods that estimate 
conductivity on the basis of electron-ion scattering cross 
sections 0, [H [H UMl ■ 

Previously 0|, the average atom Kubo-Greenwood ap¬ 
proach was systematically compared to QMD data for a 
number of elements and mixtures for a range of tempera¬ 
tures and densities. The comparison revealed reasonable 
agreement in some cases, but some serious disagreement 
in others, both qualitative and quantitative. In this pa- 
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^ Hereafter referred to as quantum molecular dynamics (QMD). 


per we adopt a similar systematic approach but using the 
Ziman-type method. While much has been published us¬ 
ing this approach, a systematic evaluation of its accuracy 
relative to QMD simulations has not, to the best of our 
knowledge, been carried out. Moreover, to generate the 
inputs (i.e. scattering potential, ionization fraction and 
ion-ion structure factor) for the Ziman-type approach we 
use the recently developed pseudo-atom molecular dy¬ 
namics (PAMD) method [Ig. PAMD has been shown 
to be of comparable accuracy to DFT-MD methods for 
equation of state, ionic diffusion coefficients and ion-ion 
structure factors [1^ [l^ , though it has limitations . 
The advantage of of PAMD over DFT-MD methods is 
that it is much less computationally expensive. As a 
result of the accuracy of PAMD, the comparison of elec¬ 
tron conductivities to QMD should be a fair test of the 
conductivity model (i.e. the Ziman-type formula), and 
not be seriously limited by inaccuracies in the underly¬ 
ing model. 


In section m we review the Ziman-Evans formula for 
electronic conductivity, and then rewrite it in terms of 
an electronic relaxation time. This allows the connec¬ 
tion and relationship of the Ziman-Evans formula to the 
recently developed expression of Johnson [2^ to be clar¬ 
ified. Moreover, we are able to generalize Johnson’s ex¬ 
pression to explicitly include the ion-ion structure fac¬ 
tor. We then give the extensions of these formulas for 
mixtures, and in the case of the Johnson formula, its ex¬ 
tension to thermal and optical conductivities. In section 
uni we summarize the pseudo-atom molecular dynamics 
method and discuss the options for coupling it to the 
conductivity formulas. Then in section IIVI the system¬ 
atic comparison to the QMD data is carried out. We also 
test the accuracy of Thomas-Fermi PAMD compared to 
Kohn-Sham PAMD for warm dense aluminum. Finally 
in section El we draw our conclusions. 
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II. EXPRESSIONS FOR ELECTRON 
TRANSPORT COEFFICIENTS 


Using this app roximation and integrating Eq. O by 
parts gives 


In this section we review the commonly used Ziman- 
Evans formula for resistivity. We rewrite it in terms of 
an electronic relaxation time in line with Ziman’s original 
derivation, and give the relationship between the direct 
approach to the conductivity (as originally proposed by 
Ziman [1^) and the inverse approach, where the inverse 
conductivity (i.e. the resistivity) is directly calculated (as 
proposed by Evans [13|). Both approaches have recently 
been used in the literature (eg. by Johnson and Faus- 
surier et al i) without discussion of their relative merits 
or relationship. Next we show how these formulations 
are easily extended to mixtures and finally, how the di¬ 
rect approach is easily extended to optical conductivity 
and thermal conductivity. We use Hartree atomic units 
in which h = rUe = e = 1. 


A. Ziman-Evans formulation 


The Ziman-Evans formula for calculating the resis- 
tivity R, which is the inverse of the conductivity, is 


OO 



0 



da{€, 9) 
d0 


Sii{q)dq 
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( 1 ) 


where e = p^/2, n* is the number density of scattering 
electrons, nj is the density of ions, Sjj{k) is the ion-ion 
static structure factor, 

=2/[I-008 0] (2) 


is the magnitude of momentum transfer squared, and 


/■ 

^ = 127r4n*r J dkk^^{k/2)V^<^<^^\kfSn{k) ( 6 ) 

where is the Fourier transform of 


B. Relaxation time formulation 


Following Ziman [2^ and Johnson [2^ we define a re¬ 
laxation time Tp 

1 1 } firr 

-=^n?- / dqq^-{p,9)Sniq) (7) 

where p = meV. This agrees with the relaxation time 
used in [2^ except that we have included a structure fac¬ 
tor, whereas Johnson’s formula has Sii{q) = 1. We can 
now write down the conductivity aoc using the Boltz¬ 
mann approach of Ziman and the resistivity R using 

M 
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As T —>• 0, df/de —)> —(5(e — ep), where ep is the Fermi 
energy. In this limit the Fermi velocity vp = (Stt^u*)^/^ 
and both these expressions give the same result 


= ^=^e 


( 10 ) 


/ = 


I 

exp(/3(e-/i*))-f 1 


(3) 


is the Fermi-Dirac occupation factor. Here /3 = 1/T is 
the inverse temperature and is the electron chemical 
potential. The differential cross section is given by 


9tT(e, 9) 
89 


y^(2Z -I- 1) sin77ie®’"P/(cos0) 
1=0 


(4) 


where pi are phase shifts and Pi{cos9) are Legen¬ 
dre polynomials. Equation o is referred to as the 
transition- or t-matrix approach to the cross section. The 
phase shifts pi in equation (j4]) are determined by solv¬ 
ing the Schroedinger equation for a scattering potential 

Yscatt^p'^ 

In the Born approximation to the differential cross sec¬ 
tion we have 


i.e. the Drude form. For T > 0 however, the results will 
differ. We refer to equation (|S]) as the direct approach 
to the conductivity, and equation ([9]) as the inverse ap¬ 
proach. It is easy to show that the inverse approach and 
the Ziman-Evans formula (equation ([IJ) are the same. 

In recovering the Drude expression (equation ((TUI) I 
from both formulations we have assumed that the chem¬ 
ical potential and electron density are related through 

nl= CTpIi/2[P9*e\ = \ j devf{e) (II) 
Jo 

where cpp is a constant [T7I| and I 1/2 is the Fermi inte¬ 
gral of order 1/2 These equations (|H]), ([U]) and m 
implicitly assume a free electron density of states. 


C. Extension to mixtures 


dcr{e,9) _ 1 yscatt(.2 
89 ^ 47r2 


To extend the Ziman-Evans formulation to mixtures 
(5) we follow ref. [s^. The extension is simply achieved by 
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the replacement: 


d(j{e, 9) 
d9 


Ns Ns 

Sii{q) Stj{q)T,ie, e)T* (e, 9) 

i 3 


( 12 ) 

where is the number species and J^i is the scattering 
amplitude for species i 


J-iU,9) = - y^(2^ + 1) sin77;e®’’* P;(cos0) (13) 
^ 1=0 


or, in the Born approximation 


J^^{e,9) 


27r 


(14) 


Xi is the number fraction of species i and Sij are the 
partial structure factors. 


D. Extensions to the direct approach 


Ziman [s^ (Ziman’s equation 9.9.7) gives the expres¬ 
sion for conductivity in terms of the integral expressions 

so that aoc = Kq. He also gives the expression for the 
thermal conductivity k in terms of these integrals 


K = 


1 

T 



{Ki? \ 
Ko ) 


(16) 


For completeness we also give the extension of equa¬ 
tion ([8]) to non-zero frequencies, but do not present any 
calculations based on it. Following Johnson |^ . 


a{uj) 


j_ r 9 ^ 3 
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(17) 


This gives the free-free contribution to the optical con¬ 
ductivity and uses a low-frequency approximation. This 
formula satisfies the conductivity sum rule 

nOO 

/ dujcj{(jj) = 

Jo 


To our knowledge, such extensions have not been de¬ 
veloped for the inverse approach, though an expression 
for the inverse thermal conductivity in the Born approx¬ 
imation has been given in ref. [32j. 




dev f(e) 


(18) 


E. Physical interpretation of the quantities in the 
Ziman-like equations 

The above Ziman-like equations (|H]) and (jH]) for resistiv¬ 
ity require input of the potential from which the electron 



FIG. 1: (Color online) An example of a calculated density of 
states (DOS) compared to a free electron DOS for aluminum 
at solid density and 1 eV. The dashed lines are the “filled” 
states, i.e. the DOS multiplied by the Fermi-Dirac occupa¬ 
tion factor. The area under either dashed curve is the number 
of ionized electrons per atom, both curves integrate to 3, as 
expected for aluminum under these conditions. The conduc¬ 
tivity formulas (equations (1S]| and (O) assume a free electron 
density of states. 


scatters the density of scatterers n*, the ion-ion 

structure factor Sii{k) and the electronic chemical poten¬ 
tial /i*. In Ziman’s original derivation [2^, jg ^ 

weak pseudopotential such that the Born approximation 
is valid and the total potential F‘°*(r) is given by the 
superposition of l/®™‘*(r) at each atomic site, i.e. 

F‘°‘(r) =^F"““(|r-i^,|) (19) 

i 

and n* is the density of valence electrons. Evans 
assumed the same form dT51) but 1/'*°“** (r) was allowed to 
be a strong scattering potential, not a pseudopotential. 
He further required that the scattering potential be of 
muffin-tin type, such that ( 7 -) = 0 for r > Rmt, 

where Rmt is some muffin-tin radius. Formally, it is 
required that no two muffin-tin potentials overlap. The 
pseudopotential in Ziman’s formulation is supposed to 
preserve the scattering properties of the full potential in 
Evans’ approach. 

For n* there are at least two reasonable choices. If the 
number of ionized electrons per ion is Z, then the density 
of these valence electron is hg = n^Z. It seems reason¬ 
able to identify n* = hg. However, the electron density 
n* is also related to the chemical potential /r*. If we take 
the limit T ^ 0 of the conductivity formulas (equations 
([5]) and ([1])), we should recover the Drude form, equa¬ 
tion m- Because of the implicit free-electron density 
of states (DOS) in equations ([U and ([9]), to recover the 
Drude form /x* must be related to n* by equation m- 
If n* = rig then we can use equation (Hip to find a corre¬ 
sponding chemical potential, /2g. However, in general the 
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FIG. 2: (Color online) An example of the potentials and 
for aluminum at solid density and 1 eV. The muffin- 
tin like extends as far as the ion-sphere radius Rws 

(~ 3 ub for aluminum at this density). The potential is 
not muffin-tin like and can have long range structure due to 
Friedel oscillations in the electron density. The inset shows 
the same potentials, but focused on the regions where the 
potentials are close to 0 . 


plasmas we consider do not have a free electron DOS. An 
example of the differences between a free electron DOS 
and a calculated DOS is given in figure [TJ The impor¬ 
tant point is that with the calculated DOS a different 
chemical potential, /ie say, is needed to recover the cor¬ 
rect number of valence electrons Z. If we use (the more 
realistic) /ig in the conductivity formulas then because of 
their implicit free electron DOS, a different electron den¬ 
sity, rig say, results from equation (ED. Thus we have 
two reasonable choices: either use hg and /ig, or rig and 
/ig. In what follows we test both choices. Finally we note 
that this ambiguity is brought about by the implicit free 
electron DOS in the conductivity formulas and would be 
resolved if a calculated DOS could be used in equations 
® and However, we do not attempt this here. 


III. UNDERLYING PHYSICAL MODEL 


sensitive to the ion-ion structure factor [s^ . Thus we 
have confidence in the inputs that it generates for the 
Ziman-like formulas. In the next section we summarize 
PAMD with the aim of describing the various options for 
how it can be coupled to the Ziman-like formulas. 


A. Summary of pseudo-atom molecular dynamics 


In PAMD the electron density of the plasma ng(r) 
is approximated by a superposition of identical pseudo¬ 
atom electron densities nf"^(r) placed at each nuclear 
site Ri 


ne(r) =Y^n^^(jr - Rij) 

i 

The total scattering potential is given by 

^ Kf^(|r - R,l) + KTK(r)] 

i 


where 



\r' — r\ 


( 20 ) 


( 21 ) 


( 22 ) 


Z is the nuclear charge of the pseudo-atom, and is 
the exchange and correlation potential. The ionic posi¬ 
tions Ri are generated in PAMD using classical molecular 
dynamics (MD), with a pseudo-atom pair interaction po¬ 
tential V{r) being calculated directly from n^^{r). Al¬ 
ternatively, if one requires only knowledge of the ion- 
ion structure factor, V{r) can be used in the so-called 
Ornstein-Zernike (OZ) equations. Solution of the OZ 
equations is much faster than using MD and is very accu¬ 
rate provided a good approximation for the bridge func¬ 
tion is used^. We have used the OZ solution for all results 
presented here. 

The calculation of V (r) has been described in detail in 
refs. [13, but it is relevant to summarize it here. One 
starts by defining an electron density ng'^’’(r), that is the 
component of (r) due solely to valence (or screening) 
electrons. In Fourier space, V(k) is given by 


U(fc) = ^Z2-nr(A:)^^ (23) 


The Ziman-like formulas require inputs of ^*'=“**( 7 -)^ 
n*, /z* and Sii{k) from an underlying physical model. 
Numerous models have been used in the literature, and 
each can be coupled to the Ziman-like formulas in differ¬ 
ent ways, depending on what one chooses for the inputs. 
Here, for our underlying model we choose the recently de¬ 
veloped Pseudo-atom molecular dynamics (PAMD) [1^ . 
This model self consistently calculates the electronic and 
ionic structure, and it agrees well with generally accepted 
DFT-MD methods on equation of state, ionic transport 
and the structure factor itself [2^,|33. Moreover, it agrees 
well with a recent X-ray scattering experiment that is 


where Z = J (PrnT'"{r) and —Cie{k)ll3 is often called a 
pseudopotential [l3j. This pseudopotential is generated 
by requiring that in the linear response regime —Cie{k )//3 
is the potential that returns the density i.e. 

nrik) = -^^Xe(fc) (24) 


^ For the bridge function have used VMHNC for single species 
plasmas [s^. Issll and HNC for mixtures [^, since VMHNC is 
not available. See f33| for further details. 
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Input choice Description 

VAA 

yPA 

•^scr 

Equation (I26II. The muffin-tin like average atom potential. 
Equation (I27II. The potential due to one pseudo-atom. 

Equation (I25II. The screened pseudopotential. 


P-e, h° 

The “physical” quantities calculated in the average atom model. The correct number of valence 
electrons is returned if fis is used with the calculated, non-free electron density of states. 

The quantities corresponding the density of valence electrons if a free-electron density of states is used. 


TABLE I: Input choices for the Ziman-like formulas that are tested here. 


where Xe is the resp onse function of the electrons (see 
equation (17) of [^). This type of pseudopotential has 
been used, along with Born approximation, to calcu¬ 
late resistivities with the Ziman-Evans formula [ij, . 
While this pseudopotential is weak by construction (due 
to the linear response relation), it is not designed to pre¬ 
serve the scattering properties of the full scattering po¬ 
tential in our case, see equation (l?7l) l. In the 

usual pseudopotential formulation (the one assumed by 
Ziman [2^) one tries to preserve the phase shifts of the 
full scattering potential, and modify the potential only 
near the nuclei. This is not guaranteed for —Cie{k)/I3 
with potentially serious consequences for the resulting 
scattering cross sections. Nevertheless, we have imple¬ 
mented this calculation and show results using it in sec¬ 
tion IIVI It is important to note that it is the screened 
pseudopotential that is used in the Ziman-like 

formulas, where 

= (25) 

e(fc)P 

where €{k) is the dielectric function dl. 

The process for calculating rig^ir) involves solving an 
average atom (AA) model. In this model, we use density 
functional theory to find the electron density (r) ^ in 
a spherically symmetric model system, in which a nucleus 
of charge Z is placed at the origin and the remaining ions 
in the system are modeled as a continuous distribution 
represented by a Heaviside step function 0(r — Rws) 
at the Wigner-Seitz radius Rws (this is called the ion- 
sphere model). A sphere of radius Rws surrounding the 
central nucleus is required to the charge neutral. The 
total average atom potential is given by 

V^^ir) = --+ f dr' - Rws)) 

r J \r-r' \ 

+H-[n{““(r)]-y-[nO] (26) 


® We retain the notation used in with the exception that 
is renamed 


where rig is the electron density in the field free region. 
The average atom chemical potential /tg is determined 
when solving the model, and is related to rig by equation 
(ED- It is a good estimation of the “physical” chemi¬ 
cal potential, as it takes into account a non-free electron 
density of states (see figure [T|). 

yAA(j,) thought of as a mufhn-tin potential, in 

the sense that Evans (30| meant. It violates the condition 
that no-two muffin-tins overlap, however this may not be 
too serious a problem since such overlaps are routinely 
tolerated when using muffin-tin potentials for X-ray spec¬ 
troscopy calculations [s^. However it is fair to point out 
that the average atom potential is not generated in the 
same way as mufhn-tin potentials are usually generated 

M- 

The screened potential (1^51) and the average atom po¬ 
tential (1261) have been widely used with the Ziman-Evans 
formula to calculate conductivities [H, [H, IH, |4l| . A 
third type of potential, that of a single pseudo-atom, has 
also been used It is important to note that while 

all these previous calculations have used potentials falling 
broadly into these three categories, differences in imple¬ 
mentation details are important and have signihcant con¬ 
sequences on conductivity calculations. In PAMD we can 
dehne a pseudo-atom potential V^^[r) as 

V^^{r) = Vri^{r) + V,^,^{r) (27) 

where the electrostatic part was dehned in equation 
(1^ . The exchange and correlation part 14^^ cannot be 
exactly separated from V).g“[ne(r)] (equation (1^1 . but a 
reasonable approximation can be written down in anal¬ 
ogy with the dehnition of the pseudo-atom electron den¬ 
sity. n^^{r) is calculated from the average atom electron 
density by removing the electron density found 

in a system with the same surrounding ions, but no cen¬ 
tral nucleus, (see [s^ for further details): 

n^^{r)=nf-^\r)-nr{r) (28) 

Therefore we define 

^ HgT[n{““(r)] - VX^\r)] (29) 

An example of V^^{r) and E'^"^(r) is given in figure [5] 
Note that does not go exactly to zero at Rws 
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Temperature [eV] 


FIG. 3: (Color online) Electrical conductivity of dense hydro¬ 
gen. All curves use PAMD-KS with a t-matrix calculation 
of the conductivity. The QMD data is from Q|. Clearly us¬ 
ing the average atom (AA) potential grossly overestimates 
the conductivity relative to QMD, while setting S{k) — 1 
demonstrates that the increase in conductivity on lowering 
the temperature is due to changes in ionic structure. 


(the mufhn-tin radius), which is consistent with typical 
muffin-tin calculations (d^ . 

In what follows we will numerically test all three po¬ 
tentials and both types of chemical potential. We will 
also demonstrate the effect of the ion-ion structure fac¬ 
tor. All of these potentials, and the structure factor are 
self-consistently calculated in PAMD. Due to agreement 
with widely accepted DFT-MD simulations on equation 
of state, ionic self-diffusion and static ionic structure, as 
well as with available experimental measurements, we 
have a high degree of confidence in the inputs to the 
Ziman-like formulas generated by PAMD. Through a sys¬ 
tematic comparison of these various inputs with QMD 
conductivity data we aim to evaluate the best choice of 
inputs, if such a choice does indeed exist. We restrict 
ourselves to comparisons where we have a high degree of 
confidence in the PAMD model. Roughly speaking this 
means temperatures > 1 eV and densities roughly equal 
to or greater than solid density. It may be possible to go 
to lower densities or temperatures but that would require 
further testing and possible improvements to the PAMD 
model. As an aid to the reader we summarize the various 
input options that are to be tested in table U 


FIG. 4: (Color online) Electrical conductivity of dense hydro¬ 
gen. All curves use PAMD-KS with chemical potential jle- 
The QMD data is from Q|. The Born calculations for both 
potentials are reasonably close to the QMD results. 


IV. COMPARISON OF RESULTS 

PAMD is a density functional theory based model. In 
it we can use either Kohn-Sham (KS) or Thomas-Fermi 
(TF) treatments of the electronic kinetic energy term. 
The KS treatment is more physically accurate and in fig¬ 
ures [3] to [10] we restrict ourselves to it. We must also 
choose an exchange and correlation potential, for which 
we have used the finite temperature local density approx¬ 
imation of [13 . Also, in figures [3l to [TOl we use the in¬ 
verse approach to conductivity only (equation (jH])). We 
address Thomas-Fermi results and the direct approach 
in sections flVBI and IIV Cl We compare to QMD calcula¬ 
tions and assume that these are close to the correct result. 
QMD is thought to be an accurate method for determin¬ 
ing the electrical conductivity in warm dense matter, but 
is very expensive computationally. There can be issues 
with convergence of these calculations [1, Ull due to the 
extreme computational expense, but in the absence of 
experiments, they are widely accepted as being the most 
accurate data available. 

In figure [3] we show the electrical conductivity of dense 
hydrogen. We have used the t-matrix method for all 
curves. Clearly, using the AA potential results in sig¬ 
nificantly too large a conductivity compared to QMD, 
though qualitatively the results with are reasonable. 
Using the pseudo-atom potential we find very good 
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FIG. 5: (Color online) Electrical conductivity of dense beryl¬ 
lium. All curves use PAMD-KS with a t-matrix calculation 
of the conductivity. The QMD data is from [3. As for hy¬ 
drogen, using the average atom potential overestimates the 
conductivity. 



FIG. 6: (Color online) Electrical conductivity of dense beryl¬ 
lium. All curves use PAMD-KS with chemical potential fie- 
The QMD data is from [3. While the Born approximation 
for the potential gives reasonable results, that for 
strongly disagrees in magnitude with the QMD data (it was 
scaled by a factor of 1/40, as indicated in the plot). 

agreement with the QMD. For we show results us¬ 
ing both chemical potential options, and find only small 
differences. However, in the physical regimes shown in 
the plot, the two estimates of n* are close. Setting the 
structure factor to be equal to 1, we see that the rise 
in conductivity at low T is no longer reproduced. In Q 
a similar effect was seen, though not explained. Those 
calculations Q used an average atom potential and no 
explicit structure factor was included, thus implicitly us¬ 
ing S{k) = 1. From figure [3] we can now understand 



FIG. 7: (Golor online) Electrical conductivity of molybdenum 
at 10 g/cm®. All curves use PAMD-KS. In the top panel 
all calculations use the t-matrix method, in the bottom all 
calculations use the chemical potential fe- The QMD data 
is from [^. In the bottom panel the I/®®’' Born results have 
been scaled by the factor indicated in the plot. 

the result of Q: The lack of the rise in conductivity at 
low T was due to the implicit S{k) = 1, while the over¬ 
estimation was because the average atom potential was 
used. 

In figured] we show the same QMD hydrogen data, but 
now test the Born approximation. Using with the 
Born approximation worsens agreement with the QMD 
data, though it is still reasonable. Using the screened 
pseudopotential 1/®'=’' with the Born approximation wors¬ 
ens the agreement still further. 

We now turn to warm dense beryllium. This is a more 
complicated atom than hydrogen and one may expect a 
comparison to QMD data to be a more robust test of 
the method. In figure [5] we show t-matrix calculations. 
Again, we see that using overestimates the QMD, 
while using with either chemical potential choice 

gives good agreement. We find that for the conditions of 
the plot, setting S{k) = 1 only has a small quantitative 
effect, though qualitatively a decrease in the conductiv¬ 
ity at low T seen in the full calculation disappears. In 
figure El we test the Born approximations. While using 
the Born approximation with results in reasonable 
agreement with the QMD, using it with 1/®'=*' leads to vast 
differences. To fit the curve on the same plot we had to 
scale it by a factor of 1/40. This result strongly suggests 














FIG. 8: (Color online) Electrical conductivity of a deu- 
terium/tritium mixture (1:1) at 10 g/cm®. All curves use 
PAMD-KS. In the top panel all calculations use the t-matrix 
method, in the bottom all calculations use the chemical po¬ 
tential fie. The QMD data is from [3. 


that, even where the Born approximation is reasonable, 
using the screened pseudopotential can lead to seriously 
inaccurate results. 

Figure [7] shows results for molybdenum. This is a much 
more complicated atom than hydrogen or beryllium and 
is a transition metal. This comparison should be an even 
more stringent test of the method. In the figure see that 
both and give reasonable agreement in mag¬ 
nitude with the QMD. Using gives poorer qualita¬ 
tive agreement at the lower temperatures, but we must 
be cautious, since it is likely that the underlying model 
(PAMD) will become less accurate at the lowest tem¬ 
peratures seen, where molybdenum approaches its solid 
phase. The two choices of chemical potential make a rel¬ 
atively large difference here, with the two curves bracket¬ 
ing the QMD data. The effect of S{k) is only important 
at the lowest temperatures. In the bottom panel of the 
figure we test the Born approximation. For both po¬ 
tentials 1/®'=’' and the results have been scaled as 
indicated, and both do a poor job of estimating the con¬ 
ductivity. Clearly the Born approximation has broken 
down here, as the large difference between the t-matrix 
and Born results for indicate. The fact that us¬ 
ing 1/®'=’’ comes closer to the data than for beryllium is 
probably due to a cancellation of errors: the Born ap¬ 
proximation underestimates croc, while 1/®'=’" leads to an 



FIG. 9: (Color online) Electrical conductivity of a car¬ 
bon/hydrogen mixture (7:9) at 10 g/cm®. All curves use 
PAMD-KS. In the top panel all calculations use the t-matrix 
method, in the bottom all calculations use the chemical po¬ 
tential fLe- The QMD data is from [^. In the bottom panel 
the Born results have been scaled by the factors indicated in 
the plot. 


overestimation. 

We now turn to mixtures: deuterium/tritium (1:1), 
figure [SI carbon/hydrogen (7:9), figure |9l beryl¬ 
lium/deuterium/tritium (1:1:1), figure Uni The results 
are broadly in line with the previous result seen in fig¬ 
ures m to [71 Using with either chemical poten¬ 

tial gives good agreement with the QMD data, with the 
choice n* = hg leading to somewhat better agreement 
in general. Using leads to an overestimation of 

the conductivity, though we find better agreement where 
the heavier element is present (the carbon/hydrogen 
mixture). In all cases the Born approximation with 
overestimates the corresponding t-matrix calcula¬ 
tion, and this overestimation becomes worse for the lower 
temperatures. Using 1/®°’’ with the Born approximation 
makes the agreement with QMD worse still, and for the 
carbon/hydrogen case it is in particularly bad agreement. 


A. Discussion of the comparison to QMD 

The comparisons in figures |3I through |Tni indicate that 
using with either chemical potential option leads 

the good agreement with QMD data. Using n* = hg 
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FIG. 10: (Color online) Electrical conductivity of a beryl- 
lium/deuterium/tritium mixture (1:1:1) at 10 g/cm®. All 
curves use PAMD-KS. In the top panel all calculations use 
the t-matrix method, in the bottom all calculations use the 
chemical potential The QMD data is from [^. 


is perhaps the slightly better option, with significantly 
better agreement being found in some cases. Using 
with the Born approximation seems to be an un¬ 
reliable method, sometimes leading to very inaccurate 
results, even where the Born approximation is reason¬ 
able. As discussed earlier, the screened pseudopotential 
is not designed for this purpose, and is not guaranteed 
to preserve the scattering properties of the full potential 
Since 1/®'=*' is generated from the screening elec¬ 
tron density n®‘^''(r), one possibility is that removing the 
oscillations in n®'^’'(r) due to core valence orthogonality, 
before generating 1/®“’’, could improve the results. We 
have tested this by using the Troullier-Martins method 
[^ . and found no significant improvement (not shown). 
Moreover, the resulting conductivity becomes strongly 
dependent on the Troullier-Martins core radius, a highly 
undesirable quality. In summary, we do not recommend 
using the U®^’' potential. 

The average atom potential leads to surprisingly poor 
agreement in magnitude for many cases. This is sur- 
jrisiM because it has been used frequently in the past 
IS. m, [13, and because it is consistent with Evans’ 
30| original formulation. The failing may be due to the 
mufHn-tin approximation. In making this approximation 
Evans’s formula is limited to non-overlapping scattering 
potentials. By forcing the potential extend only as far as 



Temperature [eV] 

FIG. 11: (Golor online) Electrical conductivity of aluminum, 
comparing Thomas-Fermi (TF) to Kohn-Sham (KS) based 
PAMD for three different densities. The PAMD-TF results 
are generally quite close to PAMD-KS, and are very compu¬ 
tationally efficient to calculate. All calculations use and 
the t-matrix approach. 

the muffin-tin radius (the Wigner-Seitz radius for 
one immediately prevents recovery of Ziman’s formula, 
which is correct in the Born approximation limit. In 
that limit, if the scattering potential is long ranged (i.e. 
longer ranged than Rws) then the muffin-tin approxima¬ 
tion will significantly change it (i.e. shorten its range). 
Only if the scattering potential is shorter ranged than 
Rws does Evans’ approach recover Ziman’s in the Born 
limit. Using is a compromise between these two ap¬ 
proaches that appears to capture the best parts of both. 
The potential is not modified by the muffin-tin approx¬ 
imation, so the Born limit can be recovered, while the 
use of the t-matrix method as proposed by Evans, allows 
strong scatterers to be present. Qualitatively, it is also 
consistent with Ziman’s original idea of having isolated 
scattering centers [^ . 


B. Thomas-Fermi results 

As noted earlier, PAMD is a density functional the¬ 
ory based model. All results presented up to this point 
have used the Kohn-Sham (KS) functional. In figure [TT] 
we compare conductivities for warm dense aluminum us¬ 
ing KS to Thomas-Eermi (TF) calculations. All results 
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Temperature [eV] 


FIG. 12: (Color online) Electrical conductivity of warm dense 
hydrogen comparing the direct approach (equation ((S])) to the 
inverse approach (equation ((U). All curves use PAMD-KS. 

use the potential and the t-matrix approach. We 

show results from both chemical potential options. Gen¬ 
erally the TF model is remarkably close to the KS based 
result. The results are closer when using rig compared 
to hg. This not surprising since the lack of atomic shell 
structure in TF leads to a continuous change in the num¬ 
ber of valence electrons with density and temperature. 
Whereas, for KS the shell structure shows up as plateaus 
in the number of valence electrons. On the other hand the 
TF estimation of the chemical potential /ie is generally 
quite close to the KS value, hence the closer agreement. 
Though TF is not as accurate a model as KS, it very 
computationally inexpensive, typically taking less than 
1 minute per density and temperature point, compared 
to roughly 30-60 minutes for Kohn-Sham calculations on 
the same modern, desktop computer. It is therefore an 
attractive option for rapidly generating tables of reason¬ 
ably accurate electrical conductivities. 


C. Direct approach and thermal conductivity 

In this section we compare the direct and inverse ap¬ 
proaches to conductivity (equations (|S]) and (0). We 
noted earlier that as T —5> 0 both approaches give the 
same result, but at high T the results will be different 
in general. In figure [1^ we compare these approaches to 
QMD data for dense hydrogen. As expected both meth¬ 



Temperature [eV] 


FIG. 13: (Golor online) Thermal conductivity and Lorentz 
number of warm dense hydrogen using the direct approach 
(equation (dl). All curves use PAMD-KS. The horizontal 
dashed lines indicate the degenerate {L = 3.29 a.u.) and 
non-degenerate {L — 1.597 a.u.) limits. 

ods give similar answers for low temperature. At higher 
T the direct approach gives much large values for the 
conductivity than the inverse approach. In comparison 
to the QMD data the inverse approach is generally closer, 
particularly at 80 g/cm^. We have also checked the com¬ 
parison to QMD for other cases (beryllium and molybde¬ 
num) but the results from the direct and inverse methods 
are similar due to the relatively low temperatures where 
QMD data is available. 

In figure [13] we show the thermal conductivity and cor¬ 
responding Lorentz number 


calculated using the direct approach (see equation (ITBl) ). 
At low temperatures the agreement with QMD is excel¬ 
lent and the degenerate limit of the Lorentz number is 
reproduced. However, for higher temperatures the agree¬ 
ment with QMD becomes poor, and the direct approach 
results in an incorrect limit of the Lorentz number. Such 
behavior was also found in Q . In that work the conduc¬ 
tivities and Lorentz number were calculated in a similar 
approach, the major differences being that the structure 
factor was implicitly set to 1 and an average atom poten¬ 
tial and chemical potential were used. We conclude that 
it is the direct approach, as opposed to the scattering po- 
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tential, that is responsible for the incorrect limit of the 
Lorentz number. Therefore the direct approach must be 
treated with caution for higher temperatures where the 
plasma is transitioning to being non-degenerate. 

V. CONCLUSIONS 

Expressions for electron transport coefficients have 
been systematically tested by comparing to accurate but 
expensive Quantum Molecular Dynamics (QMD) simu¬ 
lations. A number of implementation options have been 
tried, including choice of scattering potential, chemi¬ 
cal potential and ion-ion structure factor, all of which 
are generated by a self-consistent plasma model called 
pseudo-atom molecular dynamics [2g. For the scattering 
potential, three types of potential were tested: a muffin- 
tin like average atom potential; the potential due to one 
pseudo-atom (which is not muffin-tin like); and a pseu¬ 
dopotential method. We found that the pseudo-atom po¬ 
tential gives consistently good agreement with the QMD 
data for a range of elements (including mixtures), densi¬ 
ties and temperatures. The average atom potential gives 
good qualitative agreement but can be significantly in er¬ 
ror for magnitude. Finally, the pseudo-potential method 
was found to be unreliable, and the reasons for this are 
discussed. The effect of two reasonable choices of chemi¬ 
cal potential was tested: first, a chemical potential that 
recovers the expected number of valence electrons with a 
free electron density of states (fie)] and second, a chem¬ 
ical potential that recovers the expected number of va¬ 
lence electrons with a calculated, non-free electron den¬ 
sity of states {fJ-e)- The effect of the choice was generally 
quite modest, but using generally results in better 
agreement with the QMD data. 

Pseudo-atom molecular dynamics is a density func¬ 
tional theory based model and can therefore use dif¬ 


ferent approximations for the electronic kinetic energy 
functional. We have compared the effect of using the ac¬ 
curate Kohn-Sham approximation to the less accurate, 
but very computationally efficient, Thomas-Fermi func¬ 
tional for warm dense aluminum. In general we find the 
Thomas-Fermi results are close to those using the Kohn- 
Sham functional, and its use therefore allows very rapid 
generation of electron conductivities that are reasonably 
accurate. 

Finally, we have tested an approach in which the con¬ 
ductivity is calculated directly and compared to the in¬ 
verse approach, where the resistivity is directly calcu¬ 
lated. Formally and numerically, both give similar an¬ 
swers for degenerate plasmas, but for higher temper¬ 
atures where the plasma is transitioning to the non¬ 
degenerate regime, they give significantly different re¬ 
sults. For dense hydrogen the inverse approach was found 
to give conductivities that agree better with QMD data 
at high temperature. We used the direct approach to 
calculate the thermal conductivity and Lorenz number 
for warm dense hydrogen and found that, while good 
at low temperature in the degenerate regime, at higher 
temperatures the method does not tend to the correct 
non-degenerate limit, and significantly disagrees with 
QMD data. As a result we conclude that the direct ap¬ 
proach should be used with caution in the non-degenerate 
regime. 
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